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Abstract 

We study the hyperbolic scaling limit for a chain of N coupled anharmonic oscil- 
lators. The chain is open and with the following adiabatic boundary conditions: it is 
attached to a wall on the left and there is a force (tension) r acting on the right. In 
order to provide the system of the good ergodic properties, we perturb the Hamiltonian 
dynamics with random local exchanges of velocities between the particles, so that mo- 
mentum and energy are locally conserved. We prove that in the macroscopic limit the 
distribution of the density of particles, momentum and energy converge to the solution 
of the Euler equations, in the smooth regime of them. 

1 Introduction 

The aim of this paper is to study the hydrodynamic limit for a non-equilibrium system 
subject to an exterior time dependent force at the boundary. We consider the most 
simple mechanical model with non-linear interaction, i.e. a one dimensional chain of 
N anharmonic oscillators. The left side is attached to a wall, while on the right side is 
acting a force r (tension). For each value of r there is a family of equilibrium (Gibbs) 
measures parametrized by the temperature (and by the tension t) . It turns out for this 
specific model that these Gibbs measures can be written as a product. 

We are interested in the macroscopic non-equilibrium behaviour of this system as 
N tends to infinity, after rescaling space and time with N in the same way (hyperbolic 
scaling). We also consider situations in which the tension r depends slowly on time, 
such that it changes in the macroscopic time scale. In this way we can also take the 
system originally at equilibrium at a certain tension to and push out of equilibrium by 
changing the exterior tension. 

The goal is to prove that the 3 main conserved quantities (density, momentum and 
energy) satisfy in the limit an autonomous closed set of hyperbolic equations given by 
the Euler system. 

We approach this problem by using the relative entropy method (cf. [10 ) as already 
done in [8] for a system of interacting particles moving in K 3 (gaz dynamics). 

The relative entropy method permits, in general, to obtain such hydrodynamic limit 
if the system satisfy certain conditions: 

A) The dynamic should be ergodic in the sense that the only conserved quantities 
that survive the limit as N — >■ oo are those we are looking for the macroscopic au- 
tonomous behavior (in this case density, momentum and energy). More precisely, 
the only stationary measure for the infinite system, with finite local entropy, are 
given by the Gibbs measures. 
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B) The limit can be obtained in the smooth regime of the macroscopic equation. 

C) Microscopic currents of the conserved quantities should be bounded by the energy 
of the system. 

We do not know any deterministic hamiltonian system that satisfy condition[X] and 
this is a major challenging open problem in statistical mechanics. Stochastic pertur- 
bation of the dynamics that conserves energy and momentum can give such ergodic 
property and have been used in [8J (cf. also [7][31[2])- We use here a simpler stochas- 
tic mechanism than in [8]: at random independent exponential times we exchange the 
momentum of nearest neighbor particles, as if they were performing an elastic colli- 
sion. Under this stochastic dynamics, every stationary measure has the property to be 
excheangeable in the velocity coordinates, and this is sufficient to characterize it as a 
convex combination of Gibbs measures (cf. [2] and [I]). 

About condition [B] it is well known that nonlinear hyperbolic equations in gen- 
eral develop shocks also starting from smooth initial condition. Characterization and 
uniqueness of weak solutions in presence of shock is a challenging problem in the theory 
of hyperbolic equations. We expect shock decrease the thermodynamic entropy asso- 
ciated to the profiles of the conserved quantities (i.e. a macroscopic decrease of the 
thermodynamic entropy). The way relative entropy work in this case is in comparing 
the microscopic Gibbs entropy production (associated to the probability distribution 
of the system at a given time) with the macroscopic (thermodynamic) entropy produc- 
tion. If no shocks are present both entropy productions are small. The presence of 
the boundary force changes a bit this balance, since one should take into account the 
(macroscopic) change of entropy due to the work performed by the force. It turns out 
that the right choice of the boundary conditions in the macroscopic equation compen- 
sate this large entropy production, keeping the time derivative of the relative entropy 
small. It would be interesting to prove similar cancellation of entropy productions when 
this is caused by shocks, as it would allow to prove the hydrodynamical limit in these 
cases. Recent efforts in this direction use different methods (cf. [3]). 

About condition [c] it created a problem in [8], since in the usual gaz dynamics the 
energy current has the convecting term cubic in the velocities, while energy is quadratic. 
This was fixed in [5] by modifying the kinetic energy of the model: if the kinetic energy 
grows linearly as a function of the velocity, the energy current will grow also linearly. 
Since we work here in lagrangian coordinates, our energy current does not have the 
cubic convecting term. This allows us to work with the usual quadratic kinetic energy. 

This paper has been partially supported by Agence Nationale de la Recherche, under 
grant ANR-07-BLAN- 2-184264 (LHMSHE) and partially supported by a grant of the 
MAThematics Center Heidelberg (MATCH). 

2 The Model and the Main Theorem 

We will study a system of N coupled oscillators in one dimension. Each particle has 
the same mass that we set equal to 1. The position of atom i is denoted by q t 6 E, 
while its momentum is denoted by pi £ R. Thus the configuration space is (R x M.) N . 
We assume that an extra particle to be attached to a wall and does not move, i.e. 
(qo,Po) = (0, 0), while on particle N we apply a force r(t) depending on time. 
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Denote by q := (51, . . . , gjv) and p := (pi, . . . ,pjv)- The interaction between two 
particles i and i — 1 will be described by the potential energy — Qi-i) of an 

anharmonic spring relying the particles. We assume V to be a positive smooth function 
which for large r grows faster than linear but at most quadratic , that means that there 
exists a constant C > such that 



lim 



V{r) 



r -mo m 



lim sup V"(r) < C < 00. 



(2.1) 



(2.2) 



Energy is defined by the following Hamiltonian: 

%v(q,p): = -^Pi +^V(g l -g l _i) 



Since we focus on a nearest neighbor interaction, we may define the distance between 
particles by 

n — qt — qi-i, i = l,...,N. 
We define the energy of particle i 6 {1, . . . , JV} as 

<N := f + V(n) 



that Hn{v, p) = Xa=i e i> where r := (ri, 



The dynamics of the system is determined by the generator 

Qn ■= L T N + 7Sjv- 
Here the Liouville operator L T N is given by 

JV „ N-l 



+ (r(t)-V'(r N )) 



dp? 



(2.3) 



(2.4) 



where we used the fact that po = 0. 

The symmetric operator Sjv is the generator of the stochastic part of the dynamics 
that exchange at random time velocities of nearest neighbor particles. For any smooth 
function /, we define the operator Ti^+i by 



Tm+i = g (/( r 'P l 



-/(r, P ; 



(2-5) 
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where p l,l+ G E is defined from pel by exchanging the coordinates pj and Pj+i 

. ( Pj if j 7^ M + 1 

= < P»+i if j = i 
[ Pi if j = i + 1 

Then Sat is defined through 

JV-l 

S N f(v,p) := (/(r,p i,<+1 ) -/(r,p)) (2.6) 

i=i 

JV-l JV-l 

= -^T^ +1 /(r,p) = 2^r iii+1 /(r,p), (2.7) 



With this choice of the noise, the three balanced quantities i.e. locally conserved, are 
given by r i ,p i ,e l . 

We define £(r,p) — (r,p, —e(r,p)) G E 2 x R_, and the partition function Z on 
E 2 x R+ by 

Z{\) = Z(Ai, A 2 , A 3 ) := / e A -«('-'P) drdp. 

and the free energy function as its logarithm: 

6(A) :=logZ(A), (2.8) 

By the condition imposed on V, this function is always finite. 

For C G K 2 x E_ we define $ : I 2 x 1_ -> R by the Legendre transform of the free 
energy function 

<&(C):= sup {r? ■ C - Q(v)} ■ (2-9) 

?7eR 2 xi + 

So that the inverse is 

6(77):= sup {t}-C(r,p)-*(C)}- (2-10) 

CgR 2 xR_ 

We denote by A(u) and u(A) := (t,p,—E) the corresponding convex conjugate 
variable, that satisfy 

A = D<!>(u) and u = DQ(\), (2.11) 
where the operator D is defined by 

for any C 1 function / : E 3 -¥ R and a := (oi, 02, 03) G E 3 . 

On the one particle state space E 2 we define a family of probability measure 

vx(dr, dp) = e x<[r ' p) - ew drdp (2.13) 

Observe that 

so we can identify u = (t, p, — E) as respectively the average distance, velocity and 
(negative) energy. We also define the internal energy e = E — p 2 /2, and it is easy to 
see that A is actually a function of t, p, e. We have also the relations 

E vx (p 2 ) - p 2 = A3- 1 := /T\ P(t, e) := E vx [V(r)\ = ^ := r 
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that identify /3 _1 as temperature and r as tension. This thermodynamic terminology 
is justified by observing that, for constant r in the dynamics, and any /3 > 0, with the 
choice A = (/3t, 0, f3) the family of product measures given by: 

JV 

v$-p,a,p){<k>dp) = Y]_V(Tp,o,p)(dri,dpi), f3 G K + 

i=l 

is stationary for the dynamics. These are the grand canonical Gibbs measures at an 
average temperature pressure r and velocity 0. 

In what follows we need also Gibbs measure with average velocity different from 0, 
and we will use the following notation: 

N 

1 = 1 

where £. := (Ci,i> Ci,a. Ci,a) ■= {n,Pu-e%)- 

In a similar way we may introduce the local Gibbs measures: For any continuous 
profile u(x), x £ [0, 1], we have correspondingly a profile of parameters X(x), and we 
define the inhomogeneous product measure 

N 

"M0 : = II e A(i/JV) ' C *- e(A(i/iV)) dr i dp ij (2.14) 

i=l 

that we call Local Gibbs measures. 

We are interested in the macroscopic behavior of the interdistance, momentum 
and energy of the particles, at time Nt, as TV — > oo. Taking advantage of the one- 
dimensionality of the system, we will use lagrangian coordinates, i.e. our space variables 
will be given by the lattice coordinates {1/7V, . . . , (N — 1)/N, 1}. 

Consequently, we introduce the (time dependent) empirical measures representing 
the spatial distribution (on the interval [0, 1]) of these quantities: 

1 N ( i \ 
T)% (dx,t) := — ^6 fa;- — J &, a (Nt) dx, for a = 1,2,3. 

i=i ^ ' 

We expect the measures r]^(dx,t), a = 1,2, 3 to converge, as N — > oo, to mea- 
sures v(x,t)dx, p(x,t)dx, —E(x,i)dx being absolutely continuous with respect to the 
Lebesgue measure and with density satisfying the following system of three conservation 
laws: 

d t x - d x p = 

( to (a:) =t(a;,0), p (x) = p(x, 0), E (x) = E(x, 0) 



d t p - d x P(t, e) = 

d t E-a a ,(pp(v,t)) = o 



p(0,t) =0, P(x(l,t),e(l,t)) = T(t) 



(2.15) 

for bounded, smooth initial data to,po, Eq : [0, 1] — >■ R and the force r(t) depending on 
time t. Here we denoted by x the specific volume, p the velocity, E the total energy 
and e := E — |p 2 the internal energy. 

We need the solutions of the system (2.151 to be C 2 -solutions. To assure this, the 
following additional compatibility conditions at the space-time edges (a;, t) — (0, 0) and 
(x,t) = (1, 0) have to be satisfied: 

limpoOr) =p(0,0) =0 , Urn P(t (x), e (x)) = r(0) (2.16) 

lim jLp(xo{x), to(x)) =0 , lim ^-P(t (x), e (x)) = r'(0) (2.17) 
x— »o dx x—ti at 

^^) = ■ li S l7^^oW,eo(x)) = r"(0). (2.18) 
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A proof of this can be adapted from Chapter 3.5 and 7 of |Bj. 

For any test function J : [0, 1] — > K with compact support in (0, 1) consider the 
empirical densities 



T)%(t,J) := {rja(dx,t); J) = -^^2 J 



(2.19) 



Our goal is to show that, starting with an initial distribution such that there exist 
smooth functions to, po and Eo satisfying 



{Vi (0, J), ^(0, J), 7^(0, J)} 



J(x)xo{x)dx, / J(x)po(x)d. 



J(x)Eo(x)dx 



(2.20) 

in probability as N — > co, then at time t G [0,T] we have the same convergence of 
1)^ (t, J), a = 1, 2, 3 to the corresponding profiles x(x, t), p(x, t) and E(x, t) respectively, 
that satisfy ( |2.15[ )-( |2?l8l ). 

Here is the precise statement of our main result, where we make a stronger assump- 
tion on the initial measure: 

Theorem 2.1 (Main Theorem). For any time t £ [0, T], denote by ^ the probability 
measure on the path space C([0,T], (M. 2 ) N ) of our process with generator NQn, and 
starting from the local Gibbs measure ) corresponding to the initial profiles uq. 



Then for any smooth function J : [0, 1] 



lim v t 

N—>co 



and any S > 



J(:r)u(:r, t)dx\ > S 



(2.21) 



where u is a C 2 -solution to the system of conservation laws ( 2.15 1 — ( 2.18 1 and < T < 
t a , t a being the time at which the solution u produces the first shock. 



Remark 2.2. As our proof is based on the relative entropy method of )10j . it is only 
valid as long as the solution to (2.151 are C 2 . Since, even for smooth initial data, 
the solution will develop shocks, we are forced to restrict our derivation to a time 
< T < t a , where t s is the time when the solution to the system of conservation laws 
enters the first shock. 

Remark 2.3. A proof for the existence of smooth solutions to the initial-boundary- 
value problem (2.151 can be found in [6]. Notice that we can rewrite the pressure P as 
a function of specific volume r and entropy s: 

P(x,s) :=P(t,t). 



Then we can rewrite the initial boundary value problem (2.151, in the smooth regime, 
in terms of the unknown r, p and s(r, c) as follows: 



9 t r — d x p = 
d t p-d x P(t,s) = 
d t s = 

where we used the thermodynamic relation 



t (x) = x(x, 0), po(a;) = p(x, 0), s (x) = s(x, 0) 
p(0,t) = 0, P(x(l,t),s(l,t))= T (t) 



(2.22) 



P(t,s) = 



dc(x,s) 
dx ' 



Hence the specific entropy s does not change in time and for any x € [0, 1] is given 
through the initial data s(x, 0) := So(x). 



(i 



In the non-consevative form, equation (2.22| reads as: 



dt\p\ -A(t,p,s)d x p =0 



where the 3 x 3-matrix A is defined by 



with c : 



;(t,s) = 






1 





dp 





dp 


at 


Os 



















)-: 


— c 









0/ 



and 



s- 1 



S:=S(t,p,«) = 



c — c 
vO 



With these notations we can rewrite (|2 . 22p in the characteristic form 
S' 1 ■ d t 



c 0\ 
-lO -c ■ S _1 ■ d x 
\0 0/ 







c(d t x - cd x x) + (d t p - cd x p) + l^(d t s - eft*) = 



c 

1 dP i 



c(d t t + cflUt) - (dtp + cd x p) + ff£(ft* + cd x s) = 
d t s = 0. 



In this way we can apply the existence proof for C 2 solutions to (2.15 I for short times 
from |Bj. 



3 The Hydrodynamic Limit 
3.1 The Relative Entropy 

On the phase space (K 2 )^ we now have two time dependent families of probability 
measures. One of them is the local Gibbs measure v^r t s constructed from the solution 
of system of conservation laws ( 2.15 H 2.18 1. We denote its density by 

JV 

g N _ TT e A(«/-iv,t)-C 4 -e(A(i/jv,t)) ^ 3 
i=i 

On the other hand we have the actual distribution, whose density /t (r,p) is the 
solution of the Kolmogorov equation: 

^f(r,p) =NG* N f t N (v,p) 

(3.2) 

fo N (r,p) =ffo Ar (r,p). 

By — L T ^* + jSn we denote the adjoint operator of Qm, which can be computed as 
Lff = —L T N . 

The relative entropy of / t with respect to is defined by 

H N (t) = f f t N log^ drdp (3.3) 
J 9t 

Our main result will follow from: 



7 



Theorem 3.1 (Relative entropy). Under the same assumptions as in Theorem 2.1 
for any time t £ [0, T] 

lim ~H N (t) = 0. 

iV-KX) N 

Notice that the relative entropy H(a\/3) of a probability measure a with respect to 
a reference measure /3 and having densities / and g respectively with respect to the 
Lebesgue measure, can be rewritten as 

H(a\p) = sup tpda - log j e v dp\ (3.4) 

where the supremum is taken over all bounded functions <p. It is easy to see that 
the relative entropy has the following properties: H(a\f3) is positive, convex and lower 
semicontinuous. It follows that for any measurable function F, any positive constant 
a and some probability measures a and /3: 

E a [F] < - ]ogEp[exp(ffF)] + -H(a\/3). (3.5) 
a a 

To see how Theorem |3 . 1 1 implies the Main Theorem, please see [3[T] 



3.2 Time Evolution of the Relative Entropy 



In this Section we will prove Theorem |3.1| 

Notice that with the choice of our initial distribution the relative entropy at time 
is equal to zero: Hn(0) = 0. The strategy is to show that for some constant C 

j t H N {t) < CH N (t) + R N (t) (3.6) 

with 

lim — / R N (s)ds = 0. (3.7) 

Then it follows by Gronwall's inequality that limjv_ ) . 00 = which concludes the 

proof of Theorem |3.1| 

The following differential inequality comes from a straightforward calculation ( [TJ 

mm)- 



NGn + ^ ) log;/,' 



ft N drdp (3.8) 



Before we proceed in the proof, we have to introduce some further notations: 



For any C 1 function F := (/i, /2, /a) 



we define 



DF(r) := ((Z?/ 1 )(a),p/ 2 )(a),(D/ 3 )(a)) 5 



with D/i(a), i = l,2,3 defined by ( |2.12 l. 

Recall that u = (r, p, —E), c := E — ^p 2 and let us denote by 



J(u) :=(p,P(t,e),-pP(r,e)) 



Mr, c) Ai(r, e) 
P 'A 3 (r,e)' P A 3 (r, e ) 



the flux of (2.151, that can be rewritten as 



(3.9) 
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with the Jacobian 



DJ(u) 








1 


) 




OP 




BP 




Or 


-Pf 


9c 


v- 









With the dual relation (2.11 1, A is solution of the symmetric system 

dt[DQ(\)] — d x [DT,(X)], 

where 

E(A) = A ■ J(D9(A)). 



Equation (3.111 can be rewritten as 



Since 

it follows that 
Since 



(D 2 Q)d t \ = {D 2 T,)d x \. 

D 2 Q(X(t,x))~ 1 = (D 2 $)(i(t,x)), 
9 t A(D 2 $) = (D 2 Ti)d x \. 



(D 2 E) = (D 2 $)(DJ(u)) 
the following system of partial differential equations is satisfied: 

d t X(t,x) = (DJ) T (u)d x X(t,x). 

Let us define the microscopic fluxes: 

Ji_i,i := (-pi-i,-V'(ri),pi-iV'(ri)) i = 1, ...,N- 1, 
Jiv,jv+i := {—pN,—r(t),p N T(t)) 



Then recalling the definition of the Liouville operator given by (2.4 1 

-^AfCi ~ J»— l,i — Ji,i + 1 

Finally let us define 

v 3 :=(0, Pi ,-p 2 /2). 



Hence with the definition of the symmetric operator given by (2.61 



S N (Cj) = -2vj + v i+i + Vj_i, j = 2, . . . , N - 1 
Sn(Cn) = -v N + vn-i, Sn(Ci) = _V 1 + V 2 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



Lemma 3.2. 

2V 

NL T N \ogg? =^d x \(^,t) ■ J i - lti + NX 2 (l,t)T(t) + a N (t) (3.14) 
i=i 

where ajv(t) is smc/i t/iat 

lim — / / ajv(s)/j v dpdrds = 

iV-+oo N J J 
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Proof. 



NL T N lo ggt JV (r,p) = 7V^A(^,t)-(J l _ 1 , l ~J l , i+1 ) 



= iv53(A(-^,t)-A(5- SF i,t))-Ji-i,i-A(0,t)-J ,i + A(l,t)-J w .A 



'N N 

Taking into account the boundary conditions on A we have 

A(0,t) ■ J ,i = Aa(0, t)V'(n) = p(0,t)F'(n) = 
and 

A(l,t) ■ Jjv.at+i = -p N Xi(l,t) - A 2 (l,t)r(t) + X 3 (l,t) T (t)p N = -A 2 (M)r(*) (3.16) 



(3.15) 



because r(t)A3(l,t) = Aj.(l, t). Since A is a C -functions, we obtain 13.141 with 

N-l 



a N ( 



(*)I = SEII j *I 



»=i 

It remains to show, that liniAr-joo f* J dv^ds — 0. This is an easy consequence of 
Lemma 13.51 □ 

Lemma 3.3. 



ftiog 5 f = £(m) T (u(^,t))a,A(A t) . ( Ci-aC^.i) 



' N 



Proof. 



By (2.111, L>9 (A(^,t)) = u(i,t), and ( |3.12| the result follows. 



Lemma 3.4. Recall the definition of the symmetric operator given by (2.61. 
1 



lim . 

JV^oo N 



N S 'n log dv^ 1 ds < lim — \— 

AT— >oo 0~N 



Hisi{s)ds, 



where a is a constant independent of N with < a < i. 
Proof. 

Sivloggf 



□ 



i 1 
X(jj,t) ■ (vi_i -2v, +v i+ i) + A( — ,t) • (-Vi + v 2 ) + A(l,t) ■ (vjv-i 



»=2 
N-l 



N 



Vat) 



^(A^,*)-^' t)+A(*±i,*)).v 4 



+ (A(^ > t)-A(^ 1 t))-v 1+ (A( 



AT 



■ J V- 1 



t)-A(l,t) -vjv 
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In Lemma|3.5|we will show that the expectation of Yli ll v >ll i s uniformly bounded 
for all N and hence, since A is in C 2 , the first term vanishes in the limit as N — > oo. 

Recall that by the entropy inequality (3.5 I, for any smooth function J and any 
a > we have for k £ {1, 



1 

N 



,N}: 



PkdVs < -jrr- log 



Since this inequality is true for any a > 0, the integral on the right hand side of the 
inequality is bounded as long as a < ^ and hence the first term vanishes as N — > oo. 
The expected value of pt can be controlled in a similar way. □ 



To conclude the proofs of Lemma |3.2| and Lemma [3.4| it remains that the expected 
values of the densities are uniformly bounded: 

Lemma 3.5. 



1 



N 

i— i 

where C is a constant not depending on N . 

Proof. By the conservation of total energy and our assumption on the initial state, we 
know that 



/iV ~ iV 

^eidv?= / J2^ d "o <CN. 
i=l •> i=l 



□ 



So far we have from Lemma |3.2| |3.3| and |3.4 

d 



dt 



H N (t) < /X>A(^,i) 



(m) T (u(^,t))(C I -u(^,t) 



1 N 



-Nr(t)X 2 (l,t) + -H N (t) + R N (t) 



dvt 



(3.17) 



where_Rjv(t) is such that < |3.7[ ) holds. 
With (3.9 1 we furthermore have 



d x \(x, t) ■ J(u(:r, t))dx = 



d ( Ai(a:, t)\ 2 {x^ t) 



dx 



\3(x,t) 



dx = r(t)\ 2 (l,t). 



This means that we can replace — A r r(t)A2(l, t) by 

JV 

with an error uniformly bounded in N, and consequently from (3.171 we have 

AT 

d 



dt 



H N (t)< I *£d x \(jj,t) 



-r\T / ~ / 1 



j»-i,« - J ( u(^,t) ) - (Djy (u(if,t)) ( d - «K^,t) 



N 



* AT 



+ -H N (t) + R N (t). (3.18) 



Our next goal is to prove a weak form of local equilibrium. In view of this we 
introduce microscopic averages over blocks of size k + 1: In what follows, for any 
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vectorfield Y» := (Y M , Y 2 ,i, Y 3 ,i) : (R 2 ) 3 R 3 we denote by Y* := (Y&, Y% u Y^), 
block averages over blocks of length k + 1, where k > is independent of AT. For 
example 

Ci* = (Cf,i,<2*.i,Cs*i):=(r?,P?,-e?):= F ^ I £ C,- (3.19) 

|i-J|<i 

These blocks are microscopically large but on the macroscopic scale they are small, 
thus iV goes to infinity first and then k goes to infinity. 

For any smooth and bounded function J : [0, 1] — > R and any bounded function 

ip : R 2 -> R, 



jv «-f 
^E^)V'(n,P i ) = ^ £ £ v(r, ;P .,)-0(^). (3.20) 

»=l ,■— i,'_ii<-i 



li-i|<f 



To apply the summation by parts formula to (3.18 1, we therefore first need to do 
some cut off in order to have only bounded variables: Let d t b '■= {\ei\ < b}, and define 

Ji-!,i := 3i-i,ilc ilh and £ := Ci^c itb , 

then these functions are bounded. 

Assumptions (|2.1[| and t]2.2[ on the potential assert that by the entropy inequality 



(3.51 with reference measure du^r t j, the error we make by the replacement of Ji_i,i 
and £ 4 by J'-i^ and £j respectively is small in A^ if we can show that j^Hn(s) — > as 

For any a > small enough 



/ 



N 

I 



= i 



-E lo s / e 16 ^x(-.t) +— — 



i = l 



ATC(fe) gjy(£) 



where lini(,->oo C(6) = 0. 

Using that of A and u are in C 2 and the summation by parts formula (3.201, we 
arrive at 



i=l 



J ) -(£>J) J (u(~ *)) (C°-u(if,s) 



' N 



•AT' 



E ^ A (]v' s)>< 



|i-i|<f 



since by the cut off, IJ'-i^ — Ci I i s bounded. 
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We now introduce further block averages: 
For some small I > 0, such that I — > after iV — > oo and IN » k, we restrict the sum 
in the last expression as follows: 

N-[Nt] 



E 9 * x (k> 8 



= [Nt] 



k + : 



\l-i\<% 



i E u(^, s) -(w) t u(A, s) e fe -u(A, s) 



AT' 

+ C(fc + iV£), 



where C is a constant not depending on JV, fc and £. The error we made devided by N 
will vanish in the limit as well since I — >■ 0. 

In what follows, we denote by Af n the set of configurations in the box {i — n,...,i + 
n} and by A m the set of configurations in any box of size m + 1. 

Let < e < £ small such eN — > oo as N — > oo. To simplify the notation, we 
introduce 

A* ■= ^JV-2[JW]+fc 

._ ^2[eJV]+* 

._ A 2[ £ JV] + fc' 
._ ^N-2[N£]+4[Ne]+k 



A" 
A" 
A*" 



Performing again a summation by parts over blocks of size 2[eN] + 1, we get 
77 / / E ^ A (>7> s ) r » TTTT E J ?-i.« /fdrdpU'ds 

J ° J i=[AT<] V Kl<! / 



J0 i=[N£] 

X / 2pT + T E * (fcTT) E J '-U fs |a 4 » drdp|A»A. + c e 

JV-[JV<] t 
i=[JVf] 1/0 

y I:|I-[i\W] + [JVe]|<j y j;|j-i|<[ £ JV] 



d x X(—,s) 



— y 
fc+i ^ 

y J:|i-[M] + [JVe]|<| 



! i /•« | \ 

V ^ tl PR! ^ rflA-dsUrdplA... 
i=[i\tt] J0 1 1 ^ j;\j-i\<[sN] I 



+ Cs. 



Also here the prize is small in N because of the cut off. 
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Let denote by 

v? ,e,t (dr,dp) := f t N ' £ ' e drdp\ A * 
the measure corresponding to the density 

1 N ~ lm 1 f* 1 
ft '"' (r,p) - t J o 2 [eiV] + i E 



(3.21) 



(3.22) 



i=[iVf] 



3;b'-i|<[eJV] 



which is a function of (ri,Pi)^ {[Ne] _ 2{Ne] _ h i ... j jv-[jv£]-2[jv e )] + § }■ T h en tne last expres- 
sion reads as 



t / a 



i:|i_[]V<] + [JV £ ]|<| 

Similar we can do the summation by parts over blocks of size 2[eiV] + 1 for the term 
containing the empirical densities: 



, r t r N[i-Nl} 
ill n - , 1 



N 



E a,A(^, S )(DJ) T (a(^, S ))r l C 



i=[Nl] 



v N ' 



"• 6,fe j IVj 

»[JV£]-[JVe] aiy s «s 



t / d x \(~,t)(DJ) T U(± } t)j C{m]-[Ne]di>? 



J 6 (C*):=%*)l{ez<6;|i-Z|<^}=J(C- 



Notice also that 



The next Theorem will be proved in Section [3. 3| Known as the one-block estimate, it 
allows to replace the averages of the functions over the microscopic blocks by a function 
of the average and it is a crucial step towards the proof of the hydrodynamic limit: 

Theorem 3.6 (The one-block estimate). Let J : [0, 1] x IR + — > R 2 be a function with 
continuous first derivative. Then 



lim lim lim 

b—too k — >oo N — yoo 



where v t ' e ' is defined by (3.211 



— E 
+ 1 z - 



l;\l-[Ni]-[Ne]\<± 



[Nl]-[Ne] 



dz>f' e ' £ = (3.23) 



With this Theorem we obtain: 



Urn Hisr(t) < lim lim lim lim lim [tN / d x X{ — , s)-u> b (C [N(i _ [Ne] , u( — , s))dO t 

AT->oo b— >oo I— >0 e— >0 k — >oo JV— ioo \ J iV 1 1 1 1 iV 

-*JVy d s \(^,s) ■ ro(C[^ ] _ [JV£] u(^,s))d I >f' e: ^ 

+ J* R N ,k,eMs)ds + H ^fids^J (3 



for some cr' > and iijv it e « i> is such that 



lim lim lim lim lim 

b— >oc £— fO £— >0 k— foo iV— foe 



RN,k,e,e,b , „ 



Furthermore we defined 



24) 



Wf,(3,u):=J (3) -J (it) and ra(j, u) := (3 - u) . 
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To simplify the notation further, let fit be as follows: 

fib(3,u) := d x X ■ Wi,(3,u) - d a X ■ 1x1(3, u) 

Hence 



£>,n 6 (3,u) = (DJ b f( 3 ) ■ d x \ - d B X (3.25) 



is equal to zero if 3 is a solution of ( |3.12[ ). Consequently: 

ni,(u,u)=0, D 3 n b (u,u) = 0. 



Rewritten in terms of fit, the right hand side of (3.24| is equal to 



tN I Cl b (C[ Ne] _ [N£] ,ii(^,s))d£>?- e ' e + I R N ,k, s ,i,b{s)ds+ f Hl *Jf > ds 



Observe that the first term is equal to 

t N—[N£] 



/ / E n b (<:lu(^,s))d^ds + CsNb, 

J ° J i=[Nt] 

where C is a constant independent of N, b and e. Applying the entropy inequality (3.5 1 
on this expression, we obtain that for some a > it is bounded above by 

1 r r f N ~ m] i ) 1 r* 

- log expia ^lb(Ci, u (lTT' s ))\ du M-.t)ds + ~ H N (s)ds. (3.26) 

Hence it remains to prove, that the first term of this expression is of order o(N). 
This will be done using the following special case of Varadhan's Lemma: 

Theorem 3.7 (Varadhan's Lemma). Let v% be the product homogenuous measure with 
marginals v\ given by (2.131 and with rate function I : K 2 x R_ — > M defined by 

I(x) := $(x) — x • A + 6(A). 

Then for any bounded continuous function F onM 2 x R_ 

lira - log / e nF{i) dvl = sup{F(x) - /(x)} 

n^oo n J x 

Proof. A proof of this Theorem can be adapted from [T] [SI IS] 



□ 



In order to apply this Theorem we arrange the sum in (3.261 as sums over disjoint 
blocks and then take advantage of the fact that the local Gibbs measures are product 
measures: 

Assume without loss of generality that k + 1 devides N — 2[N£] — |, then 
N-[Ne] 



i=[N£] 



H b (Clu(^,s))= J2 E ^ 6 (C-,u(^, S )), (3.27) 



where B [Ne] := + 1) + [N£] + |; q G jo, . . . , - ^§ 2 ||. In this way, for any 

fixed j, the terms in the sum over i £ B™ depend on configurations in disjoint blocks, 
thus the random variables 

T j n i ,(C l fc ,u( — ,S)) 
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are independent under B s which is product. Consequently, if we apply the Holder 
inequality , we obtain 



~ [ log / f7 exp^cr T 3^b(Ci, u (-^, s )) \ dvx(., 3 )ds 



< -n^Tsf E log / exp L(fc + 1) £ Tj U b (Ct tt(4 a)) } 



a( ' je{-#,...,|} I i eB [Jvf] J 

1 JV-[W] ft r r . i ) 

{-|,...,|}'^' s ' 

^ ' i=[JV£] J J 1 J 

Then, since all the functions in this expression are smooth and the family of local 
Gibbs measures converges weakly, we obtain that 

i N ~ [m] /"* f f i } 

^^ (k + l)Na J£ J l0S J «p{-(fc + l)^(^,u(-, S ))jd^,, s) d S 

= ^ o J J l°g J exp|(fc + l)CTrj 6 «* ; ,u(x, S))} dvl {x . s) dxds. 



So now for each x £ [0, 1], the distribution of the particles in a box of size k is given 
by the invariant Gibbs measure with average u(x, s). such that we can apply Theorem 
|3.7| on this product measure to obtain that the last expression is equal to 

~[ f sup{an b ( 3 ,u(x,s))-I(s)}dx. (3.28) 
a Jo Jo i 

To conclude Theorem |3.1| it thus remains to show that this is equal to zero. Since I and 
fit, are both convex, since both functions and their derivatives are vanishing at 3 = u, 



it follows from assumption (2.1 1 and (2.2 1 on the potential that <rfi(3,u) < /(j) for a 



small enough. Hence there exists a a such that the last expression is equal to zero. 
This concludes the proof of Theorem |3.1| 
Since 

4;H N (t) < CH N (t) + R N ,k,e,t,b(t), 
for some uniform constant C, it follows by Gronwall inequality that 

H N (t) < H N (0)e Ct + f R N ,k,e,eAs)e C{t - s) ds 
Jo 

< e Ct (h n {0) + R N ,k,e,e,b{s)ds 
Hence the claim follows, since 

lim lim li 

b—>oo £— >0 £— »0 k—*oo N—too J n N 



inn lim Mm Mm lim / ^ N , k ^A b ^ s ) g s _ q 



3.3 The one block estimate (Theorem 3.6) 
3.3.1 Tightness 

For a fixed k> and i G {[N£], . . . , N[l - £]}, let 

*+§ 

^i e ' k (dr,dp) := f% e ' k Yl drid P' 

i=I-4~l 
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be the projection on the configurations in a block of size k around site i of the measure 
fi t > e > defined by (3.211. The density corresponding to t>tf' 1 is given by 

f t T"--=ft N ' e Xi-S-i (3-29) 

We have the following 



Lemma 3.8 (Tightness). For each k > 2 fixed , the sequence {y t f' )jv>i of probability 
measures is tight. 



Proof. We need to prove that for each n > 0, the expected value 

\i-i\<i 

for a constant C independent of TV (and n). 

Since this expectation depends on configurations only through 
we can write 

|i-I|<§ |i-!|<| 



(3.30) 
(3.31) 



In the following we define by := f^ Yii=i drdp the probability measure with density 
ft defined by the time average 



ff(dr,dp) := - / f"ds 



Using the definition (3.31 of f t ,e ' , the right hand side of (3.31 1 is equal to 

N-[N£] \ 

]V E 2[iiVTTT ^ r,jf| A ^ + J drdp 



(3.32) 




AT ^ 2feiVl + _ 



( 



j;\j—m\<sN 



\ 



N-[Ni] 



— y 

k + l 

\ !;|i-m|<5 J 

( 



f7\ 



A 2eN + k 



drdp 



N ^ / 2[eJV] + 1 . ^ + i 



m = [JVf] 



m|<eJV 



e; /t drdp 



!;|!-i|<§ 



< Cx /( 



[JV(l-!+e)]+i 



iV 



^ e m )ft N drdp 



where Cx is a constant independent of N . This inequality is true, since the last expession 

■f] 

last expression is bounded above by the total energy which again is uniformly bounded 

□ 



averages out all the e m from for m 6 {[N(£— e)] — |], . . . , [N(l— £+e)] + 1]}. Hence this 

•rgy 

(see proof of Lemma 3.5 1. This concludes the proof of (3.301. 



Lemma |3.8| asserts that for each fixed k > 2 and each fixed i there exists a limit point 
Vt'i of the sequence (v^- e ' k )N>i. On the other hand, since the sequence {v1'i)k>2 forms 
a consistent family of measures, by Kolmogorov's Theorem, for k — > oo and e — > 0, there 
exists a unique probability measure v on the configuration space {(ri,pi); 6 z £ (R 2 ) 00 }, 
such that the restriction of v on {ijj,Vi) j^^i-h i+ k] £ (R 2 ) fe+1 } is v t f. 
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3.3.2 Proof of the one-block-estimate 



Let us define the formal generator Q of the infinite dynamics by 

g-.= c + js, 

with the antisymmetric part 

d_ 

,9rj dr J+ 

and the symmetric part 



(3.33) 



(3.34) 



(3.35) 



In section \3A\ we will prove the following Proposition: 

Proposition 3.9. Any limit point v of v t l e ' satisfies the following properties: 

(i) it has finite density entropy: there exists a constant C > such that for all subsets 
ACZ 

H\A\ (m|a|<^ ,«) < C|A|, 
(li) it is translation invariant: For any local function ip and any j £ Z, 

xp d/j, = J (Tjifi) dfi 

where Tj denotes the spatial shift by j on the configurations. 

(Hi) it is stationary with respect to the operator Q: For any smooth bounded local 
function if) 

' (gip)d^ = o. 



/< 



With this Proposition, we can apply the ergodic theorem from [2J: 



Theorem 3.10 (Ergodicity). Any limit point v of i 



*(dr,dp) is a convex combina- 



tion of Gibbs i.e there exists a probability measure a(d\) on R x E + such that 
v(dr,dp)= / a(dA)JJ — !— - e XCi dr t dpi. 



The proof of Theorem 3.10 is contained in [2], see also Chapter 2 of jT] for more 
details. The idea of the proof is the following: With Proposition |3.9| one can prove that 
v is separately stationary for C and <S. This implies that the distribution of momenta 
conditioned on position v(dp\r) is exchangeable. It is here where we need the noise. 
For details of this result please see Chapter 2.3 of [T]. In Chapter 2.2 of the same notes 
[T] there can be found a detailed proof of how this implies |3.10| 



With Lemma |3.8| and Theorem |3.10| it will turn out, that |3.6| is an application of 
the law of large numbers: 



Proof of Theorem \3.6\ 

Let N* := [Nl] - [Ne]. The measure 



in expression (3.231 can be replaced 



by v t ffl since the configurations inside the integral depend on the configurations only 
through N* - | - 1, . . . , N* + f . Thus it is enough to show that for each b and each £ 



lim lim lim 

£— >0 k— too N— >oo 



E 



•1-1,1 



l;\l-N*\<i 



du^f,: k < o. 
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By Lemma 3.8 and (ii), Hi of Proposition 3.9 the left hand side can be rewritten as 



lim lim . 

£— >0 k— >oo / k 



dv 



By Theorem [3710] v is a convex combination of Gibbs measures. Taking the conditional 



expectation with respect to Q N * we are left to prove 
/ 

k + 



lim lim 

£ — >0 k — >ca It$2 



\ 



[Je AWCi dridpi a(dA) = 0. 



Because of the cut off the expression inside the integral is bounded hand hence applying 
the dominated convergence Theorem it remains to prove 



lim lim 

£ — >0 fc^OO 



-!- V j?_u-J 6 (» 
fe+i ^ ' v 



|i|<4 



[]e XWC 'W Pi = 0. 



But this is just the law of large numbers and converges to if k — > oo and holds in the 
limit as k — > oo 

□ 



3.4 Proof of Proposition 3.9 



We start with some notations: 

• By v^- n := ft'i" n!i-2!_i dridpi we denote the reduction of the measure to 
the box {i — ^ — 1, . . . ,i + §} around site i and of size n + 2. Thus its density 



ft J ls s iven b y 



/•iV.n rN i 

Jt,i '-—It l{i-f -!,...,;+§} 



• By i^' n := /jfi™ rii=i-a_i dridpi we denote the reduction of the measure i/^, 
defined in ( 3.32[ |, to a box around site i and of size n + 2 with density 



• Finally we define the density 

t,i ( r >P) : = 



2[eN] + 



T} + 1 E T il 



t,i ■ 



\j—i\<sN 



With these notations, f t ' e ' reads as: 



N-[Nl] 

ft (r,p) = 2^ fta (rfr.dp). 



With the introduction of the local averages it is easy to show that the limit point v 
are translation invariant and stationary with respect to Q: 



19 



3.4.1 Translation invariant stationary measures 
Lemma 3.11. The measure v is translation invariant. 

Proof. Let tp be a local function. Since (/ t J ^' e,fe )jv is tight, we only need to prove that 
for each z we have 

lim lim / (tp — r z ij)) f t J e ' dvdp = 

k— >oo iV— >oo J 

For a fixed i the integral is equal to: 

N-[Ni] 

/ W> - r^f^drdp = ± Y, r*m N f- eN+k (dr, dp). 

J i=[Ni] J 

Then 

fit- T^)f»r N+k drdp = — ±— J®- T.*) Y, -J^ N+k drdp 

= 2[ £ iV] + 1 / ? (Tj ^ " ^'+^)it V ! ' 2£jV+fc drdp 



j;|i-j|<eiV 



1 



2[eiV] + 1 

of— 1 



r J'^ - r ^ ] ft* drdP 

\j;\i-j\<NE j;\i-j+z\<Ne I 



converges to as Af approaches oo. Thereby u is translation invariant. □ 
Next we prove that v is invariant with respect to the formal generator Q defined by 



13.331 



Lemma 3.12. The measure v is stationary in time with respect to the generator Q — 
C + jS, that means for any bounded smooth local function ^>(r, p) 



Qi\> dv = 0. (3.36) 
Proof. Let ip be a local function. By the same arguments as in the proof of Lemma 



13.11 1 it is enough to show that 



^L^m+il £ r j{ gmT N+k drdp = 0. (3.37) 
For a fixed i define the spacial average tjj := 2 [ £ n]+i ^2j-\i-j\<Ns 7 jV , i then with 

j;|j-i|<iVe 



we can rewrite the integral of (3.371 as 



(SV0/^ 2£jV+fc drdp = i / / |c;r)/;\/iv/p./. 



tN 



s 

ii 

9s 



ds. 



We conclude the proof by observing that this expression converges to if N — > oo, 
since ip and hence tf> is a bounded function. □ 
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3.4.2 Entropy density 

Let 

V(t0,q,p) (dr, dp) := Y\ v (Tpfi,p)(dn,dpi 



and 



Then we obtain the following Lemma: 

Lemma 3.13. The limit point v has finite entropy density, that means there exists a 
constant C > such that for all subsets A k 

H(u k \u^ M) ) < C\A k \. 

In particular there exists the limit 

H{v\u T p >0 ,p) = lim —^—H(y k \v% refitP) ) = sup H(is k \v k T>0tl3) ). 



Proof. By Lemma 3.8 the sequence {p t f' )jv is tight. By the lower semicontinuity of 
the relative entropy we have 



lim H. 

iV— >oo 



■ / ~N,E,k\ k \ rr / e,k\ k \ 

k \y t ,i \ v (tP,O,0) J = H k \ Vi \V(tP,0,P) ) . 



where for each i, the limit point v\ ,k is the restriction of v to t he box Af . Consequently 
v\' — u\ A k = v k is translation invariant since by Lemma ', 



Hence 



3.11 



v is translation invariant. 



Hk (v k \v k T pfl t p) \ — H K k (y\v( T p,o,p)) ■ 
One can prove that Hn is superadditive in the following sense (see for example [T]): 

N - 2\N£\ f A jv,e,fei k \ /n qq\ 

— Hk \ V t,i P(r/3,0,« ) ■ (3-38) 



rr f~N,e,t\ N-2N* +2k+l\ 
i*N(l>t r(T/S,0,|8) ) 



On the other hand we will prove the following bound on entropy in Lemma |3.14| 



H N {p?\v? TPi0ilS) )<CN 



(3.39) 



for some constant C independent of N. By these two results we obtain 
N - 2[N£] („ N , eM k ' V V 

i=[Ni] L " J ' j;\j-i\<[eN] 

< CN 

and the Lemma follows. 



"rP,0,p))<^ 2 [eN] + 1 ^ ■ H A^f]+fc(A iV |l'(r,0 l y8)) 

AT - 2[iW] 



A 7 



□ 



To complete the proof of Lemma 3.13 it remains to show (3.391 : 
Lemma 3.14. If satisfies 

Hn (i^ I <,o,«) < CN 
for some uniform constant C > 0, then for any N G N 

H N (v? 1 1^,0,/?)) < CiV - 
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Proof. Denote by g* the density of the invariant measure t'jv^o^a) with respect to the 
Lebesgue measure. Then we obtain: 



d „ ( N , N \ f dff f t N 

H N \ Vt I ^(r,0,/3) ) = J ~Qf l°g drd P 



dt 



= J ft N NQ N ^ log dv^ PfiiP y 

= y NL T N ^du^ Tf)i0tl3) + -yN y / t JV (r,p)S , j V log/ t JV (r,p)drdp. 

Since the reference measure is invariant with respect to Qn, the first term of the last 
expression is equal to zero. Hence, with the elementary inequality 

.P 



Q log " < 2yfa(y/]l - y/a), 



we arrive at 



j t H N {v?\v( Tfii0) ) < 7 iV y 2^/f(r,p)^^//f(r,p)drdp 

= -4 7 iV Jj2 (tm+i/A^p)) drdp < 0. 

The result follows by integrating in time and by the convexity of the relative entropy. 

□ 
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